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Abstract 

The system of an atom couples to two distinct optical cavities with phase 
decoherence is studied by making use of a dynamic algebraic method. We 
adopt the concurrence to characterize the entanglement between atom and 
cavities or between two optical cavities in the presence of the phase deco- 
herence. It is found that the entanglement between atom and cavities can 
be controlled by adjusting the detuning parameter. Finally, we show that 
even if the atom is initially prepared in a maximally mixed state, it can 
also entangle the two mode cavity fields. 
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I. INTRODUCTION 



Quantum entanglement was first introduced by Einstein, Podolsky and Rosen in 
their famous paper in 1935 [1]. It lias been recently recognized that entanglement can 
be used as an important resource for quantum information processing [2]. Entangle- 
ment can exhibit the nature of a nonlocal correlation between quantum systems that 
have no classical interpretation. However, real quantTim systems will unavoidably be 
influenced by surrounding environments. The interaction between the environment 
and quantum systems of interest can lead to decoherence. It is therefore of great im- 
portance to prevent or minimize the influence of environmental noise in the practical 
realization of quantum information processing. In order to prevent the effect of de- 
coherence, several approaches have been proposed such as quantum error-correcting 
approach [3] or quantum error-avoiding approach [4,5]. 

Instead of attempting to shield the system from the environmental noise, Plenio 
and Huelge [6] use white noise to play a constructive role and generate the controllable 
entanglement by incoherent sources. Similar work on this aspect has also been consid- 
ered by other authors [7]. In this paper, we investigate an atom couples to two distinct 
optical cavities with the phase decoherence and show how the entanglement between 
atom and cavities or between two optical cavities can be generated in the presence of 
the phase decoherence. In section II, we study the system with phase decoherence by 
making use of the dynamic algebraic method [8] and find the exact solution of the mas- 
ter equation for the system. The exact solution is then used to discuss the influence of 
the phase decoherence on the probability of occupation in ground state. In section III, 
we use the concurrence to characterize the entanglement between atom and cavities or 
between two optical cavities by means of the exact solution for the system. It is shown 
that the entanglement between atom and cavities can be controlled by adjusting the 
detuning parameter. Finally, we show that even if the atom is initially prepared in a 
maximally mixed state, it can also entangle the two mode cavity fields. A conclusion 
is given in Section IV. 

II. SOLUTION OF AN ATOM COUPLES TO TWO 
DISTINCT OPTICAL CAVITIES WITH PHASE 
DECOHERENCE 

We consider the situation that an atomic system is surrounded by two distinct 
optical cavities initially prepared in the vacuum state. The Hamiltonian for the system 
can be described by [6], 

H = ^„ataW&+y(|e)(e|-|^)(^|)+^„(a|e)(^|+at|^)(e|)+^,(6|e)(^|+6t|^)(e|), (1) 

where |e) and \g) are the excited state and the ground state of the two- level atom, 
ouq is atomic transition frequency, ga(b) is the coupling constant of the atom to cavity 
modes a{b), and a (a"!"), b {b'') are the annihilation (creation) operators of a mode of 
frequency cUa and b mode of frequency cUb, respectively. In Ref.[6], Plenio and Huelge 
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use white noise as the actual driving force of the system and study numerically the 
entanglement between two optical cavities for the system in the resonant case. Here, we 
investigate analytically the entanglement between atom and cavities or between two 
optical cavities with phase decoherence by making use of the dynamical algebraical 
method. To reduce the complexity, we consider the case of uia = oji, = u. It is easy to 
verify that there exists two constants of motion in Hamiltonian (1), 

Ki = ^(t/a« a + 9b0^o) + -^{a}h + ah^) + , 

9 9 ^ 

K, = \{glb^b + gWa) - ^{a^b + ab^), (2) 
9 9 



where g = y gl + gf- It is easily proved that the operator Ki and K2 commute with 
Hamiltonian (1). We then introduce the operators as follows 

^ ^ {9aa + 9bb)\e){g\ ^ _ {gaa^ + gbb^)\g){e\ 



9VK1 9VK1 

^o = ^(|e)(e|-|^)(^|). (3) 

It can be shown that the operators Si {i — 0, ±) satisfy the following commutation 
relations 

[5o,^±] = ±^±, [5+,5-]=25o, (4) 

where 5*0 and S± are the generators of the SU(2) algebra [9]. In terms of the SU(2) 
generators, we can rewrite Hamiltonian (1) as 

H = u;{K, + K2-^) + ASo + ^/^(^+ + S.), (5) 

where A = a;o — a; denotes detuning. In this paper, we consider the pure phase 
decoherence mechanism only. In this situation, the master equation governing the 
time evolution for the system under the Markovian approximation is given by [10] 

f^--^[H,p]-l[H,[H,p]], (6) 

where 7 is the phase decoherence rate. Noted that the equation with the similar form 
has been proposed to describing the intrinsic decoherence [11]. The formal solution of 
the master equation (6) can be expressed as follows [8] , 

p{t) = Y.^^M\t)p{m'%ti (7) 

fc=0 ^■ 

where p(0) is the density operators of the initial atom-field system and M^{t) is defined 

by 

M\t) = H''exp{-iHt) exp{-^H^). (8) 
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By means of the SU(2) dynamical algebraic structure, we obtain the explicit expression 
for the operator M'^ 

M\t) = l[/+(i^i,/r2)]'^expH/+(/ri,/r2)t]exp[-^[/+(/ri,i^2)]'] 

+ \[f-{K,, K2)f expH/_(i^i, K2)t] exp[-^[/_(i^i, K2)f] 

- [f4K^, K,)]' expH/_(Ki, K2)t] exp[-^[/_(Ki, K^)]']}, (9) 

where 

/±(Xi, K2) = cuiKi + - i) ± ^Q(Xi), n{K,) = (A^ + 4g^K^y/^, 

H^nt = ga{a\e){g\ + a^g){e\) + g,{h\e){9\ + &^l^?)(e|). (10) 

We assume that the cavity fields are prepared initially in vacuum state |00), and the 
atom is prepared in the excited state |e). The time evolution of pit) can be written as, 

Pit) = + ^ + (1 - ^) cosQt exp(-|Q2)]|00)(00| ® |e)(e| 

+■^{^[1 -cosQtexp(-^Q^)] +isinQtexp(-^Q2)}|00)((^| ® \e){g\ 
+^{^[1 -cosQiexp(-^n2-)] _isinniexp(-^n2-)j.|(^^^QQ| ^g, |^^^g|^ 

\ L \ L Z Z 

+ ^[1 - cosf)texp(-|f]2)]|(^)((^| ® \g){g\ (11) 

where 

|(^) = -(^„|10)+^,|01)), n = (A^ + 4/)V2. (12) 

The |</7) in Eq.(12) is a single-photon entangled state. Recently, much attention has 
been paid to investigate the preparation of the single-photon maximally entangled 
state [12]. It is noted that when the two coupling coefficients ga = gb, the state \ip) 
is nothing but a single-photon maximally entangled state. We then show that if a 
projective measurement on the atom in the {|e), \g)} basis is made, the atom will be 
projected on the ground state \g) with the probability Pg in the case of A = 0, 

= ^[1 - cos(2^t) exp(-27/i)]. (13) 

If the measurement result is \g). the two distinct cavity fields arc in the single-photon 
maximally entangled state ^(|10) + |01)). In Fig.l, we plot the probability Pg as the 
function of time t for different values of phase decoherence rate 7. It is shown that 
if the decoherence rate 7 is zero, the two distinct cavity fields are in the maximally 
entangled single-photon state at the time t — with unit probability. 
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III. THE ENTANGLEMENT BETWEEN ATOM 



AND CAVITIES OR TWO OPTICAL CAVITIES 

In order to quantify the degree of entanglement, several measures [13] of entangle- 
ment have been introduced for both pure and mixed quantum states. In this section, 
we adopt the concurrence to calculate the entanglement between atom and cavities or 
between two optical cavities with the phase decoherence. The concurrence related to 
the density operator p of a mixed state is defined by [14] 

C{p) = max{Ai - A2 - A3 - A4, 0}, (14) 

where the Aj {i = 1, 2, 3, 4) are the square roots of the eigenvalues in decreasing order 
of magnitude of the " spin- flipped" density operator R 

R = p{ay® ay)p*{ay® ay), (15) 

where the asterisk indicates complex conjugation. The concurrence varies from C — Q 
for an unentangled state to C = 1 for a maximally entangled state. 

We flrst investigate the quantum correlation between the atom and cavity modes. 
If we deal with the two cavity modes as system B, and the atom as system A, then 
pit) in Eq.(ll) can be thought of as the density operator of a two-qubit mixed state. 
In the basis |11), = |00) ® |e), |10), = |00) ® \g), |01), = \v) ® |e), |00), = \v) ® \g), 
the explicit expression of the concurrence Cab describing the entanglement between 
the system A and system B can be found to be. 

Cab = ^{^[1 -cosniexp(-^02-)|2 ^g.^2j^^g^p^_^^j^2^|i/2_ ^^g^ 

Prom Eq.(16), we can see that the detuning A plays a key role in the quantum correla- 
tion between the atom and cavity modes. If the decoherence rate 7 is not equal to zero, 
the concurrence Cab remains in the value 2g\l^\/Q? in the limit t 00. In the strong 
coupling case, i.e., Qa^Qb ^ A, the concurrence Cab of the stationary state p(oo) is ap- 
proximately \I\\/(2g). On the other hand, in the large detuning limit, the concurrence 
Cab of the stationary state is approximately 2g/\A\. In Fig. 2, the concurrence Cab is 
plotted as a function of the time t and decoherence rate 7. We show the concurrence 
Cab as a function of the detuning parameter A and the decoherence rate 7 at a fixed 
time in Fig.3. In the limit t 00, concurrence Cab is plotted as a function of the 
detuning parameter A in Fig. 4. From Fig.4, we can see that Cab increases with the 
detuning A parameter. This means that the entanglement between the atom and the 
cavity fields can be controlled by adjusting the detuning parameter. Now, we turn our 
discussion to the resonant case, i.e. A = 0. In this case. Cab = \ sin(25ft)| exp{—2g'^^t). 

In Ref.[15], it has been proved that for any pure states of three qubits 1, 2 and 3, 
the entanglement is distributed following the inequality for the squared concurrence 

^12 + ^13 — ^1(23)' i^'^) 

where Ci,(23) is the single-qubit concurrence defined as the concurrence between the 
qubit 1 and the rest of qubits (2,3). For any mixed states of three qubits 1, 2, and 3, 
there is analogous inequality for the squared concurrence as follows 

C!2 + Cls<{C')T(^), (18) 
5 



where (C^)^23) is the minimum of average over all possible pure state decomposition 
of the three qubits mixed state [15]. 

In the present paper, we may expect that the pair entanglement between the atom 
and the a{b) mode cavity field is determined by the coupling coefficient Qaigb)- It is 
easy to prove that there exist the simple relations, 

Cls-Cl + Cl (19) 

and Ca/Ci, = Qa/ 9b) where Ca(Cb) is the concurrence describing entanglement between 
the atom and the a{h) mode cavity field. Thus, our result is in agreement with that 
obtained in Ref.[15]. 

Next, we investigate the entanglement between light fields of two distinct cavities. 
By tracing out the degree of freedom of the atom in density matrix p{t) in Eq.(ll), we 
obtain the reduced density matrix pBif) describing the two light fields as follows, 

P^^^^ = + §^ + " cosOtexp(-|i7^)]|00)(00|, 

+ ^[1 - cosl]texp(-^Q2)]|^)(^|. (20) 

Then, the concurrence Cb characterizing the entanglement of two light fields can be 
derived as 

Cb = (21) 

From Eq.(21), we can see that the concurrence Cb is equal to zero at time t — 
2mi/Vt,{n = 0,1,2...) in the case of 7 = 0. At these specific time, the two cav- 
ity modes have no pair entanglement. However, in the case with 7 7^ 0, the two 
cavity modes is always entangled for the time t > 0. In Fig.5, we plot the concur- 
rence Cb as the function of time t and damping rate 7. Prom Fig.5, we see that the 
entanglement between the two distinct light fields increases with the phase decoher- 
ence rate 7 within the time range 2nn/\/A'^ + Ag'^ < t < (2n + |)7r/v^A^ + Ag'^ or 
(2n + l)^l/^/WTW <t<{2n + 2)7r/VA^ + V {n = 0, 1, 2, ...). The concurrence 
is displayed as a function of the phase decoherence rate 7 for three different values 
of the detuning parameters at fixed time in Fig. 6. The stationary state entanglement 
of the two cavity modes measured by concurrence is 4gagb/ (^'^ + 4:g'^) . This means that 
the stationary state entanglement achieves its maximal value 1/2 in the resonant case 
with ga = Qb- 

Finally, we discuss how much entanglement between the two mode cavity fields 
can be achieved if the initial atom is prepared in a thermal state and the cavity fields 
are prepared in the vacuum states. We assume that the initial atom is in the state 
Pa(0) = 5\g){g\ + (1 — (5)|e)(e|, where < 5 < 1, and the cavity fields are still in the 
vacuum state |00). Our calculation shows that = (1 ~ 5)Cab and C'^ = (1 — 5)Cb- 
This means that even if the initial atom is prepared in a maximally mixed state 
\\9){9\ + |k)(e|, it can still entangle the two mode cavity fields. In this case, the 
concurrence equals | in the steady state for A = and Qa — 9b- 
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IV. CONCLUSION 



In this paper, we investigate analytically the entanglement between atom and 
cavities or between two optical cavities with phase decoherence by making use of the 
dynamic algebraic method. It is found that the entanglement between atom and cav- 
ities can be controlled by adjusting the detuning parameter. Finally, we show that 
even if the atom is initially in a maximally mixed state, it can also entangle two mode 
cavity fields initially prepared in vacuum state. The approach adopted here can be 
employed to investigate the entanglement between two optical cavities mediated by a 
two-level atom in those cases, in which the two mode cavity fields are initially prepared 
in another separable states. 

ACNOWLEDGMENT 

This project was supported by the National Natural Science Foundation of China 
(Project NO. 10174066). 



7 



References 

[1] A. Einstein, B. Podolsky and N. Rosen, Phys. Rev. 47, 777(1935). 

[2] D. P. DiVincenzo, Science 270, 255 (1995); L. K. Grover, Phys. Rev. Lett. 79(2), 
325 (1997); J.I. Cirac and P. Zoller, Nature 404, 579 (2000); M. A. Nielsen and I. 
L. Chuang, Quantum Computation and Quantum Information (Cambridge Uni- 
versity Press, Cambridge, 2000). 

[3] P. W. Shor, Phys. Rev. A 52, 2493(1995); A. M. Steane, Phys. Rev. Lett. 77, 793 
(1996); E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 

[4] L.-M. Duan and G.-C. Guo, Phys. Rev. Lett. 79, 1953 (1997); P. Zanardi and M. 
Rasetti, Phys. Rev. Lett. 79, 3306 (1997); D. A. Lidar, L L. Chuang, and K. B. 
Whaley, Phys. Rev. Lett. 81, 2594 (1998). 

[5] H. M. Wiseman and G. J. Milburn, Phys. Rev. Lett. 70, 548 (1993); D. Vitah 

[6] M. B. Plenio, S. F. Huelga, Phys. Rev. Lett. 88, 197901 (2002). 

[7] A. Beige, S. Bose, D. Braun, S.F. Huelga, P.L. Knight, M.B. Plenio, and V. Vedral, 
J. Mod. Opt. 47, 2583 (2000); M. S. Kim, J. Lee, D. Ahn, P. L. Knight, Phys. 
Rev. A 65, 040101 (R) (2002). 

[8] Jing-Bo Xu and Xu-Bo Zou, Phys. Rev. A 60, 4743 (1999). 

[9] Y. Wu, Phys. Rev. A 54, 4534 (1996); Jing-Bo Xu, Xu-Bo Zou and Ji-Hua Yu, 
Eur. Phys. J. D 10, 295 (2000); Y. Wu and X. Yang, Phys. Rev. A 63, 043816 
(2001). 

[10] C. W. Gardiner, Quantum Noise (Springer- Vcrlag, Berlin, 1991); W. H. Louisell, 
Quantum Statistics Properties of Radiation (Wiley, New York, 1973). 

[11] G. M. Milburn, Phys. Rev. A 44, 5401 (1991) 

[12] M. D. Lukin, A. Imamoglu, Phys. Rev. Lett. 84, 1419 (2000); J. C. Howell, J. A. 
Yeazell, Phys. Rev. Lett. 85, 198 (2000); A. Rauschenbeutel et al, Phys. Rev. A 
64, 050301 (2001); 

[13] W. K. Wootters, Quantum Inf. Comput. 1, 27 (2001), and references therein. 

[14] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[15] V. Coffman, J. Kundu, W. Wootters, Phys. Rev. A 61, 052306 (2000). 



8 



Figure Caption 



FIG.l. The probability Pg as a function of the time t for various values phase damping 
rate : 7 = 1 (Solid), 7 = (dash), 7 = 0.01 (dot) and 7 = 0.05 (dash dot) with 
i/a = fifft = 1 and A = 0. 

FIG. 2. The concurrence Cab as a function of the time t and the phase damping rate 
7 for gr^ = = 1 and A = 5. 

FIG. 3 The concurrence Cab as a function of the detuning parameter A and the phase 
damping rate 7 for ga — Qb — ^ and t — IQ. 

FIG. 4 The concurrence Cab of the steady state as a function of the detuning param- 
eter A for Qa — gb — i and 7 = 0.1. 

FIG. 5 The concurrence Cb as a function of the time t and the phase damping rate 7 
for g^^l and A = 0. 

FIG. 6 The concurrence Cb as a function of the phase damping rate 7 for various 
values of the detuning parameter: A = (Sohd), A = 1 (Dash) and A = 2 (Dot) 
with t — 2 and ga — Qb — 1- 
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